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Abstract 
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the known determinant formula of norms and scalar products. 
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1 Introduction 

The doping of spin chains has directly observable consequences in experiments. For example 
the magnetic susceptibility of a spin- 1/2 Heisenberg chain made of Sr2CuO^ shows a strong 
dependency upon the oxygen content [2]. In one spatial dimension the impurities cut the 
chain and serve as effective boundary fields. This breaks the translational invariance and 
makes the local magnetization dependent on the lattice site. However, by introducing an 
additional reflection algebra Cherednik [4] and Sklyanin [16] showed the open Heisenberg 
chain to remain integrable. For diagonal boundaries, 
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+ \a\ coth C + a z L coth £ + ] sh rj , (1.1) 



it can even be solved by the coordinate [1] or the algebraic Bethe ansatz [16]. L is the num- 
ber of lattice sites and aj, a = x,y,z are the usual Pauli matrices. The boundary fields are 
parametrized by the complex numbers £ and r] is the crossing parameter of the model enter- 
ing the anisotropy A = chrj. The distribution of Bethe roots in the ground state depends on 
the boundary fields. For example anisotropics | A| < 1 allow for at most two purely imaginary 
Bethe numbers besides real roots [17]. From the coordinate Bethe ansatz such states con- 
taining imaginary rapidities are usually termed 'boundary bound states'. This terminology 
is related to the exponential decay of phase factors. The ground state for anisotropics A > 1 
can be found in [9]. 

To calculate the S^-magnetization in the ground state as an expectation value in the 
framework of the algebraic Bethe ansatz one could make use of the inverse problem Wang [20] 
solved in terms of a mixture of the reflection and the Yang-Baxter algebra. However, this 
method makes use of the translation operator for which the Bethe states are no longer 
eigenstates. This difficulty Kitanine et al. overcame by reducing the problem to the algebra 
of the periodic chain where the inverse problem [12] is only expressed in terms of this algebra. 
So they determined the action of local operators on Bethe vectors in the representation of 
the reflection algebra and thus were able to apply the trigonometric generalization [10] of 
the rational determinant formula [21] for scalar products. Additionally they succeeded in 
simplifying the combinatorial part of the local ^^-magnetization by introducing a generating 
function [11]. The result was a multiple integral representation for its state average value. 
The integral representations is linked to expressions from the vertex operator approach [8] 
and was derived for the ground state of (11.11) described by Bethe root densities which are 
valid in the thermodynamic limit L — > oo of the half infinite chain. 

Treating finite chains arises from the subtle observation that on a formal level integrable 
systems of finite length share much of the properties of systems at finite temperature: In 
the first case the logarithmic derivative of the eigenvalue of the usual row-to-row transfer 
matrix [3] determines the ground state energy of the system of length L. In the second 
case the free energy in the thermodynamic limit is connected to the leading eigenvalue of 
the quantum transfer matrix [18, 19] at temperature T. As recently shown by Damerau et 
al. [5], utilizing this analogy the method of non-linear integral equations [13,14] describing 
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the quantum transfer matrix could easily be applied to the usual row-to-row transfer matrix 
of the finite-length periodic Heisenberg chain to calculate high-precision data of correlation 
functions. Extending this approach to the reflection algebra of the open XXZ spin chain will 
be the scope of this article. 

The paper is organized as follows. In section [2] we will review the reflection algebra. Then 
in section [3] we give a non-linear integral equation for an auxiliary function accounting for 
the lowest lying state of f 1 1.1 1) with zero magnetization. In section H] we show our main result, 
how to combine the known determinant formula with the auxiliary function to represent 
normalized scalar products in terms of multiple integrals. Section [5] is devoted to an example. 



2 Integrable Boundary Conditions 



Sklyanin's construction [16] of integrable systems involving boundaries is valid for a gen- 
eral class of integrable systems characterized by an i?-matrix of difference form R(X, /i) = 
R(X — /x)e End(V (g) V) (V is a vector space with dim V G N ) which not only satisfies the 
Yang-Baxter equation 



R 12 (X - //) Ri 3 (X - v) R 2 s((J> 



R23 (a* - v ) #13 (A - v) R 12 {X - fi) 



(2.1) 



but also several conditions such as symmetry with respect to the permutation operator P on 

V <g> V (Px ®y = y®x), 

R{\) = PR{X)P , (2.2) 
unitarity involving some complex function p(X), 

R(X)R(-X) = p(X) (2.3) 

and crossing unitarity for another complex function p(A), 

R l1 {X)R h (-A - 277) = p(A) . (2.4) 

The parameter r\ characterizes the i?-matrix and the superscript tj denotes the transposition 
with respect to the jth space in the tensor product V ® V . Here we will need the well-known 
6- vertex model solution 



R{X) 



/sh(A + r?) \ 

shA shrj 

sh i] sh A 

\ sh(A + 77)/ 



(2.5) 



of the Yang-Baxter equation (12. ip . It generates the Hamiltonian of the spin-| XXZ chain 
with 

p(A) = sh(A + 77) sh(-A + r/) , p(A) = sh(-A) sh(A + 2r/) . (2.6) 

Each solution R(X) of the Yang-Baxter equation fixes the structure constants of a Yang- 
Baxter algebra 

it! 12 (A - /i)7i(A)T 2 (/j) = T 2 {p)T 1 {X)R l2 {X - p) (2.7) 
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with generators T a p (X), a, (3 = 1,2; where T\(A) = T(A) ® J 2 , T 2 (A) = 7 2 <g> T(A) are the 
embeddings of the monodromy matrix T(A) with 2x2 unity I 2 . 

Sklyanin's construction of open spin chains is based on the representations of two algebras 
and defined by the relations 

R 12 (X - /i)T 1 ( - ) (A)i? 12 (A + m)^2 H (^) = Tt\v)Ri2(\ + ^)T 1 H (X)R 12 (X - /i) (2.8) 



i? 12 (-A + / x)T 1 (+)<1 (A)i? 12 (-A -fi- 2r ? )T 2 (+) ' 2 ( /U ) = 

= r 2 (+)42 ( / i) J R 12 (-A-/i-2r ? )r i (+) ' 1 (A) J R 12 (-A + yU ) . (2.9) 

We shall call and T^' right and left reflection algebras respectively. The transfer matrix 

t(X) =trr (+) (A)T ( - } (A) (2.10) 

as a trace in auxiliary space is the central object under consideration because it generates 
with [t(X), t(fj,)] = a commuting family of operators. The explicit construction of integrable 
open boundary conditions for models arising from the Yang-Baxter algebra starts with the 
2x2 matrix 

^« = ih^hA( Sh(A + -sh(A-o) =/2 + ,T ' tanhACOth5 (2 ' U) 

originally found by Cherednik [4] and cast into the form (12.111) by de Vega et al. [6]. It 
constitutes the c-number representations K^^X) = \K(X + £ + ) and K^~'(X) = K(\,£~) 
of the reflection algebras with the obvious properties 

tr#(A,0 = 2, K ( -\0)=I 2 , trfsT (+) (0) = 1. (2.12) 

The Hamiltonian (11.11) involving the interaction of two neighbouring sites is connected, 
up to a factor, to the first derivative of t(X) by looking at the expansion t(X) = 1 + 2XH + . . . 
around the point A = 0. Considering local L-matrices building up the two representations 
rW(A) = L L (X) ■ ■ -L M+1 (X) and T<->(A) = L M (X) ■ • -L X (A) of then by construction 

TH(A) = T^(X)K^(X)T^\-X) 
T^ t (X) = T^ t (X)K^ t (X)(T^- 1 )\-X) 

are representations of the reflection algebras such that the normalized transfer matrix 

t(X) =tiK^\X)T(X)K^(X)T- 1 (-X) , t(0) = 1 (2.14) 

is independent of the factorization of T(A) = T^ + \X)T^ (A) . Thus we are free to choose 

T (+) '(A) = T t (X)K ( - +)t (X)(T- 1 )\-X) , T { -\X) = K { ~\X) . (2.15) 



5 



In order to gain more symmetric arguments and to avoid inconvenient scalar functions after 
applying the inversion formula 

r ''^ P^(Rf ' () -»' (2 ' 16) 

it is instructive to define the new object U ( - + \\ + r]/2) := T^(X) (d q T)(—X — r]/2) consisting 
of 

E/(+)*(A) = T\\-r 1 /2)K^ t (\-T 1 /2)a y T(-\-T 1 /2)a y . (2.17) 
It is still a representation of the right reflection algebra with a 2 x 2 matrix in auxiliary space, 

( >A) 0(A), 

\c(\) p(a); 

The quantum determinant (d q T)(X) is the central element (Casimir) of the Yang-Baxter 
algebra (12.71) . With the one-dimensional projector P{ 2 onto the antisymmetric (singlet) state 
in the tensor product V ® V of auxiliary spaces the definition reads 

(d q T)(X) = tr 12 PrMX - r]/2)T 2 (X + 77/2) 

= A(X + 77/2) D(A - 77/2) - 5(A + 77/2) C(X - 77/2) . "" 

Here, the trace tri2 is to be taken in both auxiliary spaces 1 and 2 of the tensor product 
V ® V and the monodromy matrix T enters with the usual representation 

T(A) = l'- 4 ^ B(A)> ] . (2.20) 
1 C(A) £>(*),' 

As an example, implying (12.1 7p the operator /3(A) can be reduced to the operators of the 
periodic chain via 

B ' A » = 2ch(A + U)shC 3h i 2 (2l)" ) [ Sh ' A " ^ + « + » B < A " "/ 2 » B '" A " ^ 

+ sh(X + r]/2-^ + )B(-X-r]/2)D(X-r]/2)] . (2.21) 

Remark. The related transfer matrix 

r(A) =tiU i+ \X)K^\X-r]/2) = t(X - -q/2){d q T){-X) (2.22) 

to the monodromy U^ + \X) reduces at the point X = rj/2 to (d q T)(—r)/2). As the Hamiltonian 
is proportional to t'(0) it is now related to the logarithmic derivative of r(A) at X = rj/2 
because o/t'(0) = A lnr(77/2) - d x \n(d q T)(-T]/2). 
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3 Auxiliary Function 

The open Heisenberg chain (11. ip is related to the fundamental representation 

T(A) = R 0L (\ -s l )... R 02 (X - s 2 )R i(\ - Sl ) 



(3.1) 



of the ordinary monodromy matrix in auxiliary space for L lattice sites each equipped with 
an inhomogeneity Sj e C. Choosing the vacuum |0) = (1)® L for the algebraic Bethe ansatz 
the expectation values of the operators A(X) and D(X) read 



t(A) = JJsh(A-s/ + 7/) , d(A) = JJsh(A 



(3.2) 



l=i 



and T(A) acts as an upper triangular matrix. Because of (12.171) the monodromy matrix 
lf( + \\) of the right reflection algebra acts as an upper triangular matrix with respect to 
|0) as well enabling the algebraic Bethe ansatz [16] with diagonal boundary fields The 
quantum determinant takes the explicit form 



(d q T)(X) = a(X + r)/2)d(X-r]/2) 



(3.3) 



and for the new transfer matrix r(A) = tr U^ + ' (X)K^ \X — rj/2) the corresponding Bethe 
ansatz equations [16] for M < L/2 Bethe numbers are given by 



sh(A i - + 77/2) sh(Aj - r + 77/2) 
sh(A i + e + -V2)sh(A i + e--V2) 



n 



sh(Aj - 77/g + si) sh(A i - 77/2 - si) 
\ sh(Xj + rj/2 + si) sh(A j + r?/2 - si) 



M 

n 

1=1 



sh(Aj - A; - r;) sh(Aj + A; - 77) 
sh(Aj - Xi + 77) sh(Aj + Xi + 77) 



(3.4) 



for all j 



A(z) 



+ 



M rendering the eigenvalue 



-1) 



[z + 77/2) sh(2^ + 77) sh(z + £ + - 77/2) sh(z + f~ - rj/2) q(z - 77) 



2ch(0 + 77/2)ch(z-r//2) sh(2z) sh£+ sh£" q(z) 

(-1) L (f)(z - 77/2) sh(2^ - 77) sh(z - + v/2) s\i(z -f + v/%) ?0 + v) 



2ch(z + rj/2)ch(z-Tj/2) sh(2z) 



sh£- 



q(z) 



(3.5) 



analytic at the Bethe roots Xj defining the function q(z) := [Yli=i s ^ 1 ( z ~ ^1) sn ( z + 
The shorthand <f>(z) : = [Il^=i s ^( z ~ s i) sn ( z + s 0] accounts for the pairwise distinct lattice 
inhomogeneities Sj regularizing combinatorial expressions in the forthcoming sections. 

Remark. Note that the Hamiltonian (II. ip corresponds to the homogeneous case Sj — > such 
that we have to perform this limit (including d\ln(d q T)(—rj/2) = 0) if we want to compare 
e.g. with results from the exact diagonalization. 
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Figure 1: Solution of 1 + a(z) = in the rational limit for the ground state of L = 4 
lattice sites: A typical distribution of (two) Bethe roots • and (five) hole-type solutions ▲ 
for positive boundary fields £ + = l.li, £~ = .9i (left panel) and a negative boundary field 
£~ = — .3i along with £ + = 4i (right panel) in the sector of M = 2. 



Let us restrict the anisotropy ch?7 of the zz-interaction to the massless case, rj = ij, 
and choose < 7 < 7r/2 for the region next to the isotropic point. By selecting the low- 
est lying state of zero magnetization, not necessarily the ground state, from Bethe vectors 
B(Xi) . . . £>(Aa/)|0) with M = L/2 some simplifications occur for the auxiliary function 



sh(z - £+ + vy/2) sh(z - f~ + ry/2) sh(2z - ry) 
sh{z + - ry/2) sh{z + £~ - ry/2) sh(2z + ry) 



r L 



X 



-j-j- sh(z - ij/2 + si) sh(z - vy/2 - s t ) 



1=1 



sh(z + ij/2 + si) sh(z + ij/2 — sj) J q(z — vy) 



q{z + h) (3.6) 



L 1% 

associated with the unique solution {Aj}^ =: {A}, 1 + a(Aj) = 0. Obviously L has to be 
even and both eigenvalue and auxiliary function are periodic in m. Because of periodicity 
the boundary parameters £ can be restricted to the complex interval (— ivr/2, i7r/2] for an 

______ L/2 

hermitian Hamiltonian fll.ip . Once a set of Bethe numbers {A} = {Xi}^ is fixed satisfying 
1 + a(Aj) = for all j = 1, . . . , L/2 there are additional hole-type solutions {x} = {Xk\k=i 
to the same equation, 1 + a(xk) = 0. These are also zeroes of the eigenvalue (13. 5p . 

A( , _ (-1) L 4>(z + 77/2) q(z - rj) sh(2z + rj) sh(z + - rj/2) sh(z + £~ - 77/2) 1 + a(z) 
^ ~ 2ch(z + r}/2)ch(z-r]/2) sh(2z) sh^+ sh£~ q(z) 

The number of holes follow from the transformation w := e 2z of 1 + a(z) into a rational 
function of w where the nominator is a polynomial of degree 3L + 4: Due to the symmetry 
a(— z) = l/&(z) all zeroes Aj and Xk appear twice with different signs and thus they are 
symmetrically distributed with respect to the origin as shown in figure [TJ Additionally the 
equation 1 + a(z) = has two trivial solutions z = and z = in/2 fixing the number of 
hole- type solutions to be L + 1. 



Compared to the case of the half infinite chain where the problem can be treated by 
root densities in the thermodynamic limit we want to pursue another way [5, 13, 14] valid 
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Table 1: Ten possible combinations of the boundary fields £ 



for a finite number of lattice sites. It turns out that the meromorphic function a(z) is suf- 
ficiently well determined by the gross properties of {A} and {x} depending on the value of 
both boundary fields. As (13. 6p is symmetric in the parameters £ ± , one has to distinguish 
between ten main cases, c.f. tabled], for the pole structure in view of £ ± , the positions of 
Bethe numbers and hole-type solutions. Indeed, the last case X is sensitive to the values of 
£^ but inverting all parameters £ — > — £ =b formally reverses the ^-direction and maps the 
region X to the cases / to W. 



By means of table [T] lets consider some examples by looking at the first column. From 
numerical evidence one observes for opposite boundary fields (region VII) L/2 real Bethe 
roots inside the strip | Im z\ < 7/2 and the hole-type solutions to lie outside of it. Additionally 
one hole-type solution seems to stick to the pole z = rj/2 — £ _ of the boundary field (figure 
[TJ right panel). This observation along with the eigenvalue and the known asymptotics is 
enough to derive a set of equations relating the second logarithmic derivatives of a and (1 + a) 
to each other determining a(z) uniquely by means of the integral Fourier transform. As this 
technique is explained in detail in [14] we may leave with the homogeneous Sj — > result 



In a(z) = 4:7] + In 



2r] + In 



'sh(2,2 — rj) sh(z — i])' 
sh(2z + rj) sh(z + rj) . 
sh(z- (£+-?7/2))sh(z 



(r-va)) 



+ 2r]L + 2L In 



valid for the region | Im z\ < 7/2 



Lsh(z+(£+-77/2))sh(z+(£- 



~sh(z — rj/2)" 




Lsh(2 + 77/2)J 


-L 



-r}/2)). 
sh(2^) ln(l + a{u)) 



(3.7) 



27ri sh(z — u + rf) sh.[z — lu — rj) 
s. The factor £ < 1 ensures the Fourier integrals to 



converge and serves in the inhomogeneous case as a convenient restriction | Imsj| < e. The 
canonical contour C is depicted in figure [2] and extends to infinity. Compared to the periodic 
chain [5] only the first and second lines of driving terms were added. 

For boundary fields exceeding i7 (region I) a hole-type solution on the real line appears 
besides the outermost Bethe root (figure [U left panel). This is due to a change in the 
description of the ground state encoded by L/2 real Bethe roots. All other L hole-type 
solutions remain outside the strip |Im^| < 7/2. The hole-type solution on the real line 
corresponds to the term 



+ 477 + In 



sh(z + x — rj) S M"2 — X ~ v) 
sh(z + x + v) S M<2 — X + V) 



(3. 
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Re 



Figure 2: The Bethe roots are enclosed by the canonical contour C for the massless case in 
its parametrization rj = ij, < 7 < tt/2. Note the symmetry a(—z) = l/a(z) mirroring 
all solutions of the equation 1 + a(z) = at the origin. The poles at u — tj/2 — £ of the 
function (1 + 0) have to lie outside the contour. 



which has to be added on the RHS of f ]3 . Tl) due to complex analysis imposing the additional 
constraint 1 + a(x) = on the auxiliary function. This is similar to the case of excited states 
in the periodic XXZ chain [15]. 

For boundary fields 7 — e > > (regions II, IV) the pole at z = r]/2 — £~ has 
to remain outside the contour guaranteed by a deformation. Applying the residue theorem 
yields the additional driving term 



2r] + In 



shO + 77/2 + t 



sh(z-3r?/2 + £- 



(3.9) 



which, besides (13.81) . has to be added on the RHS of (13. 7p for that case. This is the only 
modification compared to case I because the structure of the root distribution of Xj and Xk 
with respect to the strip | Imz| < 7/2 is unchanged. Nevertheless approaching with the pole 
z = 7]/2 — C,~ zero from below the Bethe root closest to the origin moves towards zero along 
the real axis (region II) . Passing the origin the pole picks up this Bethe root and pulls it (up 
to exponential corrections with respect to the chain length) along the positive imaginary axis 
until the upper part of the canonical contour is reached. Because of this imaginary Bethe 
root the corresponding states are termed 'boundary bound states' (region IV). 

Especially in the XXX limit when the parallely oriented boundary fields (region VI, two 
imaginary Bethe roots) become strong enough to significantly arrest the outermost spins of 
the chain the considered 'boundary bound state' with M = L/2 Bethe roots and a total 
magnetization of zero refers no longer to the ground state. Here the almost fixed boundary 
spins can be regarded as effective boundary fields for a spin chain with two lattice sites and 
one Bethe root less. For the XXZ chain this effect already sets in for the regions V, VI, VIII, 
IX but depends on the values of the anisotropy ch 77 compared to the boundary fields. 

Numerics suggest for all cases I to IX with zero magnetization to have L/2 Bethe numbers 
within the contour C as in the considered examples above. Thus all possible forms of driving 
terms with respect to hole-type solutions and poles of the boundary fields are given. The 
non-linear integral equation for the auxiliary function can then be fixed if one considers 
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I II III IV V VI VII VIII IX 

hole-type solution • • • • • • 

left boundary pole (— ) • • • • • 

right boundary pole (+) • • • • • 

Table 2: Driving terms for the possible combinations of boundary fields £ 



table [2] as a building block. Here • marks the additional driving terms ( 13.81) and ( 13.91) which 
have to be added on the RHS of ( 13. 7ft for each single case. This accounts for the hole-type 
solution x inside the canonical contour imposing 1 + a(x) = and the boundary fields £ ± . 

In the following two sections we shall derive our main result valid for distributions of 
Bethe numbers in the strip | lmz\ < 7/2 according to the left panel of figure [H For this 
reason, we have to introduce a closed contour C similar to C but excluding all hole-type 
solutions, especially the one on the real line closing the set of Bethe numbers. 



4 Integral Representation for the Determinant Formula 



To calculate scalar-valued expectation values of local operators a nice combinatorial result 
for the Bethe-eigenvectors YlbLi S(Aj>) 1 0) of the open XXZ chain applies. The key element is 
the inversion formula 



m— 1 



[ II (Msi) + D( 8j ))]T^8 n )[]l (A( Sj ) + D^))- 1 ] (4.1) 



for the standard basis (e. 



S^Sp, at site m from [12]. Because it is written in terms 
of the monodromy T(A) its action on a Bethe state with (I2.2ip can be computed by YBA. 
Then some Bethe numbers Xj are replaced [10] by pairwise distinct lattice inhomogeneities 
(k = v/2 + s k to regularize the expressions 



M 



(4.2) 



6=1 



Here we have {//&} = {Ajl^UlCfelfcLi and the summation is taken over certain subsets a m of 
{1, 2, . . . , M + m}. For a local operator at site m only the first m inhomogeneities s\, . . . , s m 
enter and their shift of 77/2, £fc = Sk + i]/2, is due to the explicit decomposition^ of B(X) in 



1 Note that the operator B(X) here, (|2.21[) . and the corresponding expression in Kitanine et al. differs by 
an overall prefactor and a shift of rj/2 in the periodic chain operators. Looking up [10] we find 



B 



Kit 



= (-1) 



L sh(2A + 7 ? ) 



sh(2A) 



sh(A - 77/2 + £+)B(\)D{-\) + sh(A + t]/2 - Z + )B(-\)D(\) 



To make use of the normalized scalar product formula (|4.3p one should always bear in mind, that the right 
Bethe vector containing B(fJ>k) gets its arguments from commutations starting with B(Xj) such that the 
prefactors in front of the square brackets cancel due to the normalization. 
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terms of the periodic chain operators. The coefficients C am can be computed generically and 
for an illustrating example to this formula see ( 15. 31) . 

Proposition [10]. For a set of pairwise distinct numbers {[ik}ff=i an d Bethe roots {A/}^ 



solving the Bethe ansatz equations (3.4 ) the normalized determinant formula for scalar prod- 
ucts including members of the right reflection algebra reads 



(oi[n a A iig(Aa)][n£ig(^)]io) 

<0|[n^iC(A a )][n^iB(A 6 )]|0> 



sh(A afe ) sh(A 



■abj 



n s h(^ a6 )sh(/7„,,)_ 



M 

n 



sh(2//j + rf) sh(2A, 



with the entry 



_Alsh(2A, + 77)sh(2^)_ 



det [H(\j,fi k 



j,k=l,...,M 



.M 



(4.3) 



(4.4) 



of the determinant, the shorthands X a b 
roots included in the functions 



sh(A j - fjt k ) sh(Xj + fjt k ) 
= X a — Xt, X a b := A a + Xb and the set {A} of Bethe 



y(z, {A}) 



sh(2; — Xj — rf) sh.{z + Xj — rf) 
y{z, {A}) := a(z - rf/2)d{-z - rf/2)s\i{z + - rf/2) sh(z + £~ - v/2) 

M 

x J Yl S H Z - Xi-7]) sh(z + Xl-T]) 
1=1 



(4.5) 



(4.6) 



Here a(A) = Yl l=1 sh(A — Si+r]) and d(X) = n^i s ^(A — sf) are the vacuum expectation values 
of the operators A(X) and D(X) of the periodic chain with inhomogeneities Sj approaching 
zero in the homogeneous limit. 



The Bethe ansatz equations (13. 4p follow from 



a[z 



sh(2z + rf) sh(z + Xj — rf) sh.{z — Xj — rf) 
sh(2,2 — rf) sh(z + Xj + rf) sh(;z — A,- + rf) 



(4.7) 



and can be rewritten Vj{Xj) = yj(-Xj), j = 1, . . . ,M allowing to recast the entries of the 
determinant in the form 



H(Xj, /ifc) 



sh(2/i fe - rf) 



Vj(fik) sh(/i fc + Xj - rf) sh(/i fc - Xj - rf) 
sh(2/i fc - rf) sh(Aj - fi k ) sh(A j + fi k ) { sh(/i fc + Xj - rf) sh(fi k - Xj - rf) 

sh(2/i fc + rf) 



sh(// fc + Xj + rf) sh(/i fc - Xj + rf) 



• (4- 
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Considering the limit \i k — > A& in the above expression to get in contact with the matrix 
elements to be calculated yields 



lim 



1 



sh(2/i fc - 77) 



x k sh(Xj + /i fc ) sh(Aj - fi k ) { sh(/i fe + Xj - 77) sh(/i fc - Xj - 77) 



sh(2/i fc + 77) 



sh(/i fc + Xj + 77) sh(/j, k - Xj + 77) 



shr7sh(2Aj) 



d In a 

iK v (Xj + X k ) - iK v (Xj - A fc ) - 8 3 k —^-(Xj 



(4.9) 



where one separately has to treat the case \i k — > Aj accounting for the Kronecker 5^ by virtue 
of l'Hospitals rule. K v is the kernel from the auxiliary function, 

Sll(277) 



KvW i sh(A + 77) sh(A - 77) ' 
Obviously all normalized expectation values (14.31) contain the elementary ratio 
det [i/>{K,Hb)] ab=1 



det [4>(Xj, A/ 



j,k=l,...,M 



det 



(Xj,X k )ijj(Xk,Hi) 



j,l=l,...,M 



(4.10) 



(4.11) 



where on the RHS _1 (Aj, X k ) denote the entries of the inverse matrix and summation over 
k is understood. Now for reshaping the RHS we closely follow [7] as similar was done for the 
periodic XXZ chain. Here the entries are 

iK v (Xj + A*) - iK v (Xj - A fc ) - 5i^(Xj) 



cj)(Xj, A fc ) 



ip{Xj,fi k ) 



sh(2A i ) 



(4.12) 



sh 77 sh(2/ifc — 77) 



sh(Xj - /ifc) sh(/ifc - Xj - 77) sh(/ifc + Xj - 77) sh(Aj + fi k ) 

sh rj sh(2/i fc + 77) 



- a(fx k )- 



(4.13) 



sh(Xj - /ifc) sh(/ifc - Xj + 77) sh(/ifc + Xj + 77) sh(Xj + fi k ) 

defining the new matrix J(Xj,fii) := </> _1 (Aj, X k )tp(X k , fii). Expressing the determinant (14.111) 
in terms of a density function is presented in appendix El and summarized in the following 

Lemma 1. For simplicity assume to be a copy of the Bethe numbers where the first n 

roots Ai, . . . , X n are replaced by some c-numbers Vx, . . . , v n thought of lattice inhomogeneities 
(j = 77/2 + Sj from the strip | Im(£j — 77/2)! < e. Then the considered determinant reduces to 



det [ip(X a , n b )] 



a,b=l,...,M 



det \d)(X», X 



J- 



det [j(A is i4)] =det 



j,k=l,...,M 



a'(Xj) lj,l=l,...,n 



(4.14) 



where ' denotes a derivative and G(X, v) is the solution to the linear integral equation 

sh 77 sh 77 



G(X,u) 



sh(A + u) sh(A + v — 77) sh(A — v) sh(A — v + 77) 

du sh(27/) 



+ 



G{u,u) 



c , 2m sh(A — lu + rj) sh(A — to — rj) 1 + a(uj) 



(4.15) 
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on the contour C . Here we already made use of = and mind the simple zeroes 

G(0,v) = G(X,7]/2) = 0. The density function shows the symmetry G(—X,v) = —G(\,v) 
with respect to the first argument A whereas the second argument v here is restricted to the 
strip | Im(z/ — rj/2)\ < e outside C . The contour C excludes the hole-type solutions \ and 
depends on the parameter e as shown in figure 0. 



To generalize the result let us introduce the disjoint union of the sets {A} = {A + } U {A~} 
and {/i} = {/i + } U {A~} and denote the cardinality of the partitions {A 1 * 1 } by |A ± |. Then 
along with the slightly modified function [10] 



s ct ({a+},iV}|{a^}) 



tt £0*+, {A}) sh(2/x+ + r,) sh(2A+) sh(2a+A+ - rj) 
l\ sh(2//+) sh(2//+ - 77) y(a+\+,{\}) S h(2\+ + v) 



x 



n 

a<b 



sh(A+ - A+) sh(A+ + A+) 
sh(/i+ - sh(/x+ + //+) 



n M—n 



nn 

L a=l b= 



sh(A+-A fe -)sh(A+ + A fe -) 
\ sh(/x+ - A^) sh(//+ + A fe ") 



the normalized scalar product 

(oii[nf = ig(A a )][n & =igM]io) 
(oi[n^iC(A a )][n^B(A b )]io) 

effectively reduces with |A + | 



^({A + },{/i + }|{A-})det 



a'(Ai 



j,k=l,...,n 



(4.16) 



(4.17) 



- = n to an n x n matrix. The set {a} with Oj = ±1 
accounts for the symmetry of the Bethe roots which can be seen from the Bethe ansatz 
equations in the form y(Xj, {A}) sh(— 2Aj — rf) — y(—Xj, {A}) sh(2Aj — rj) for j = 1, . . . , M 
leaving S a unchanged. 



5 Generating Function of the Magnetization 

For an illustrating example we shall now apply the integral representation of the scalar 
product formula (14.171) to a generating function of the S^-magnetization. Note that we will 
assume the case of one hole- type solution on the real line accounting for < ^/i < n/2. 

Proposition [11]. Corresponding to one of the simplest non-trivial one-point functions in 
the open spin chain is the one-parameter generating function 

m m 

Qmiv) = [n + e ^-))] [n {A( Sj ) + Disj))- 1 ] (5.i) 

3=1 3=1 

of the longitudinal magnetization 

1-a: 



2 f =D m ^(Q m (y,))|^ =0 . (5.2) 
It includes a discrete derivative D m u m = u m — w m _i on the lattice and a continuous one with 
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respect to ip. Its action on a Bethe state reads 

M 



Q m (v)[U B ^]\°) =H E E [nE] det [ M K A M* 

1=1 n=0 |A+|=n |C+|=n j=l<Tj=±l 



j,fe=l,...,n 



X 



X 



n n U( + \+\ 1 



M-n 



5; 1 ({A + },{C + }|{A-})[lI B (Ca)] [ II W)]|0> ( 5 - 3 ) 

a=l 6=1 



wii/i i/ie known matrix 
M(\ j: n k ) 



shrj 



sh(Xj - n k ) sh(A j - Hk-rj) 
e^ sh.n 



+ 



+ 



sh(\j - pLk) sh(\j - pLk + v) 



sh(Xj - A+ - rj) sh(A 3 - - + 77) 

jLi sh ( A i - ^ + - ^) sh ( A i - K + ^) 



n 



(5.4) 



from the generating function of the zz- correlation and the function 



W-({u},{z}) 
W({u},{z}) 



y sh(w, + e--77/2) 



n 



ria,6=l Sh (^ + W « - ^) 



n a<6 sh(2; a + 2 fe - 77) sh(u; a + cj 6 - 77) 



sh(^ - uj a -i]) sh(z b -u) a + rj) 
chti sll (^ a -Ub-i]) sh(z a - z b + 77) 



n 



(5.5) 



(5.6) 



picking out the left boundary with £ to start counting the lattice sites. All inhomogeneities 
entering the generating function are included within the expressions [7] 



b(A) 



sh(A - Q) 
Ai- sh(A - + 77) 



n 



(5.7) 



The function S a already appeared in the scalar product formula, whereas Uj = ±1 accounts 
for the symmetry of the Bethe roots. 

Here the combinatorial part is expressed by the set of all ordered pairs ({A + }, {A~}) of 
fixed cardinality |A + | and \( + \ respectively indexing the sums. Switching to the normalized 
scalar product the expectation value of the generating function can be written in terms of 
the density function G(o~jXj, v) = <jjG(Xj, u), 



(<2 m (y?)) = 



(oii[n^iC(A a )]g m (y)[n£ig(A b )]io> 
(oii[nf=iC(Aa)][n^B(A b )]|o) 



ee e [nE] n^: 



n=0 \\+\= n \C+\=n j=l <7j=±l 



t\ b '(Ca + ) sh(2C+-77) 



(5.8) 



xW_({a + X + },{( + })det M^+A+C) 



det 



j,k=l,...,n 



a'(At) 



j,k=l,...,n 
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The last step now is to get rid of the explicit dependence on Bethe roots by integrals 
according to the following 

Lemma 2. Let f{oj\, . . . ,u n ) be a complex function, symmetric in its arguments and equal 
to zero if any two of its arguments agree up to a sign. Furthermore if it is analytic on and 
inside the simple n-fold contour {C') n and shows a simple zero at ujj = to compensate the 
first order pole of the auxiliary function 1/(1 + a) then 



E E ••• E 

|A+|=n a +=±l 



±1 



nr=i"'W 



nl 



n 

i=i 



, 2ni 1 + a(uji) 



f{u u ...,u n ) . (5.9) 



For lattice inhomogeneities from the strip | Im(^ — T]/2)\ < e all poles with respect to the 
variable Uj of the density function G(ujj,(k) he outside the contour C. The singularity of 
the function W_ ({<*;}, {(}) at u>j = rj/2 — £~ is balanced by the simple zero of 1/(1 + fx(uj)) 
and because of the density function G(0, v) = the expression in (15. 8p meets along with the 
determinant property the conditions of lemma [2J 

However, the same technique can be applied to the C + -summation with inhomogeneities 
(k in the vicinity of rj/2 and thus outside C. For a function / with the properties from above 
except the simple zero at ujj = the corresponding integrals read [7] 



E 

15+1=' 



/(Cl~; • • • ? Cn) 

nr=i 



i 

n! 



n 

i=i 



dzi 1 
r 2nib(zi) 



f(zi, . . . , z n ) 



(5.10) 



With respect to the integrand the contour V lies in the strip | Im(z — r}/2)\ < e and 
surrounds all inhomogeneities Ci> • • • > Cm- 111 addition due to the simple zero G(X,r]/2) = 
the point rj/2 can be enclosed enabling the homogeneous limit (k — y v/2 yielding 

Proposition. Letm be a site index counted from the left boundary and consider the functions 
b(A), W-({u}, {z}), M{ujj,z k ) andG{ujj,z k ) according to (E2D, (l5T5]l . (E3D and (HT5]) . Then 
the multiple integral representation of the generating function reads 



m 1 r- n p 

n=0 v > L 1=1 JL 

x det [M( 



dzi 1 



' c , 2-k'i 1 + a{uj) J T 2-k'i b(zi 

G(ujj, z k 



WL({w},{«}) 



det 



J J,A:=1, 



sh(2z fc - r/) 



j,k=l, 



(5.11) 



T/ie contours for the massless case are depicted in figure^ and in the homogeneous limit the 
auxiliary function a(z) is determined from the non-linear integral equation (13.7)) . 



Thermodynamic Limit 

Rewriting the density function (14. 15)) as the sum G(X, v) = G + (X,u) — G + (X,r] — v) the 
partial density G + satisfies for | Im(z/' — rj/2)\ < e 

n+(x n_ shr/ . f sh{2r]) G + (u,u' 

Lr {A,V ) — --T7T ,, A -TT\ ~, , „\ + 



sh(A — i/') sh(A — z/ + r/) y C / 27ri sh(A — to + 77) sh(A — lu — rj) 1 + o(u;) 

' (5.12) 
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Figure 3: In the massless case rj = ij, < 7 < n/2 the contour C is limited by the hole-type 
solution x as depicted in the left panel for a small £<1. The lattice inhomogeneities (k he 
in the vicinity of rj/2 outside C and are counterclockwisely surrounded by V (right panel). 



In the limit of infinitely many lattice sites L — > 00, thus x ~ > 00 an d C — > C, the auxiliary 
function a is dominating for 7 < the upper part of the contour C, 

sh(A - i7/2) 



lna(A) ~ 2L- 



(5.13) 



sh(A + ij/2) ' 

such that only the lower part remains (figure HI left panel). Clearly, the validity range of 
the variable v' extends according to the pole structure of the driving term and the density 
G + (X, v') — > i7rp(A, u') satisfies in this limit 



p(AX) + 



i sh(2r])p(u, v') 



sh 7] 



2tt sh(A — to + rf) sh(A — lu — rj) n sh(A — u') sh(A 



v) 



(5.14) 



Remember, the variable v' takes the values of v and rj — u where v is located in the vicinity 
of rj/2. Applying the thermodynamic limit to the generating function one is directly led to 
the result of Kitanine et al. [11] where 7 < is assumed, 

dzi b(ui) 



<«-(»>» = Ess? n 



n=0 



1=1 



Cr 



x det M(a;j, Zk) 



r 2ti\ b(zi 
det 



W-.({u,},{z}) 

p(u j: Z k ) - p(Uj,Tj- Z k ) 

2sh(2z fc -rj) 



- j,k=l,...,n 



(5.15) 



j,k=l,...,n 

Cd consists of the real line and an additional counterclockwisely closed contour around the 
pole u = rj/ 2 — £~ dependent on the value of the boundary parameter £~ . The condition 
— I7I/2 < Im£~ < for this additional contribution to be included can easily be seen from 
figure IH right panel. 



6 Conclusion 



In this paper we stated for the open XXZ spin chain a non-linear integral equation for an 
auxiliary function accounting for the ground state of the model in regions /, 71, III, IV 
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Im 



m 



Im 



ilyl 
2 




ilyl 
2 



Re 



Figure 4: In the thermodynamic limit only the lower part of the contour C remains (left 
panel). Moving it towards the real axis the poles at ujj = rj/2 — £~ of the function W- must 
not be crossed (right panel). 



and VII. In the other cases a neighbouring state of zero magnetization is described which 
depends on the value of the anisotropy and boundary fields and reduces just to region VI in 
the XXX limit. Our formulae are derived for the Bethe ansatz solvable case of S^-conserving 
boundaries and are valid for a finite number of even lattice sites. In the main part we showed 
how to combine this method with the scalar product formula of Bethe vectors. This yielded 
a linear integral equation whose solution builds on the well known determinant representing 
scalar products. As an example we derived for a certain generating function of the S z - 
magnetization a multiple integral representation showing the correct thermodynamic limit 
and matching the result from exact diagonalization for small lattice sites. 

For the derivation we assumed a distribution of Bethe solutions having one hole-type 
solution on the real line. Thus we had to choose a closed contour C for the integral repre- 
sentations of the determinant formula and generating function differing from the canonical 
contour C of the auxiliary function. But this is already the general case. In the simplest 
issue of field parameters > ^/i > — tt/2 and n/2 > £ T /i > (regions VII, VIII and 
IX) no 'holes' have to be taken into account such that the canonical contour even applies 
for the integral representations. Therefore we expect this to be a good starting point for 
numerical considerations. Interesting as well would be a refinement of the regions V, VI, 
VIII and IX for the Bethe ansatz with reference to the alternation in the description between 
the ground state and its zero-magnetization neighboured state which both degenerate in the 
thermodynamic limit. 

In this article we restricted our derivation to the massless case rj = i'y with < 7 < tt/2. 
For the region n/2 < 7 < 7r the imaginary parts of the contours take values ±(7 — 7r/2 — e), 
e <C 1 to ensure the hole-type solutions to be located outside the contour due to i7r-periodicity. 
We would like to note that the results can be extended to the massive regime of real crossing 
parameters rj > 0, where the contours have to be redefined according to figure [5) 
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Forschungsgemeinschaft under grant numbers SE 1742/1-2 and FR 737/6 . 
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Figure 5: In the massive case 77 > the canonical contour C with a small e <C 1 is depicted 
in the left panel. The lattice inhomogeneities (k he in the vicinity of rj/2 outside C and are 
counterclockwisely surrounded by the closed contour T (right panel). 



A The Density Function 



In order not to overload the notation lets consider the set {ni}f=i = {vj)™ =l U {Xi}f£ n+1 . 
Then the expression under the determinant of the RHS of (14. lip reads in components 



^ Aj ' ^ = sh(2A.) [ K v( X J + X k) ~ K v( X j ~ X k) J(\k, vi) - 



J(Xj, vj) d In a 



sh(2Aj) dz 

solved for the part containing the logarithmic derivative of the auxiliary function 

At 



(Ai) , (A.l) 



d In a 

J ( A i> ^)^-(Ai) = iJ2 \ K n( X i + X k) ~ k v( X j - x k) J( x k, 



dz 



k=l 



+ 



shrj 



sh rj 



sh(ui - Xj) sh(^ - Xj - rj) sh(z/, + Xj) sh(z/; + Xj - rj) 

sh 77 sh 77 

sh(z// — Xj) sh(vi — Xj + rj) sh(z// + Xj) sh(z/ ; + Xj + rj). 



(A.2) 



Strikingly the definition of J(A, v) according to (1A.2j) is compatible with the constraint 
J(X,v) = —J(—X,u) and the properties a(z) = l/a(—z), a'(Xj) = a'(-Xj) of the auxiliary 
function. Considering 

F{Xj : vt) := J{Xj : v{) d -^{Xj) (A.3) 

for arbitrary arguments Xj the analytic properties ResA=±^ -F(A, z/) = 1 + a(v) beside the 
single zero F(0, v) = are known from the RHS of (1A.2|) such that 

du sh(2r/) F(u, v) 



c , 2ni sh(A — lu + rj) sh(A — lu — rj) 1 + a(a>) 

sh?7 



F(X,v) 



+ (! + *(")) 



+ 



sh?7 



.sh(A + v) sh(A + v — rj) sh(A — v) sh(A — v — rj) 



(A.4) 
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holds where all Bethe roots Xj and the variable v except the hole-type solution \ on the real 
line are supposed to lie inside C . Redefining F(X, v) := —(1 + a(u))G(X, v) we are led to 

. sli77 sh?7 
G{X,v) = —— — — L - + 



sh(A — v) sh(A — v — rj) sh(A + v) sh(A + v — rj) 

+ f du MM G(u,u) ^ 

Jc 27ri sh(A — u + rj) sh(A — c<j — rj) 1 + a(c<j) 

Unfortunately in our case v should be a lattice inhomogeneity Q taken from the strip 
| Im(z/ — r)/2)\ < e. Thus G(X, v) is due to an additional residue and a(£j) = the solution 
to the linear integral equation (I4.15p . The considered determinant (14. lip is calculated from 



det^fA^Afe)!., , J ' L -'^ ^ 1 - 1 ; 



i,fc=i,...,M 



a (ui))G (Xj, vi) 



j,l=l,...,n 



and making again use of Cl(£j) = for the variable V\ — Q — r]/2 + Sj reduces the expression 
( 1A.6I) to ( 14. 14f) . Clearly explicitly using a(Q) — changes G(X, v) for the argument v away 
from (j compared to the original definition. But because G(X, v) is only used in combination 
with simple poles and the residue theorem (c.f. ( 15. 10D ) just its unchanged value G(X,Q) a t 
the lattice inhomogeneity Q is relevant. 
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